Effect of Dipolar Interaction in Molecular Crystals 
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We investigate in this paper the ground state and the nature of the transition from an orientational 
ordered phase at low temperature to the disordered state at high temperature in a molecular crystal. 
Our model is a Potts model which takes into account the exchange interaction J between nearest- 
neighbor molecules and a dipolar interaction between molecular axes in three dimensions. The 
dipolar interaction is characterized by two parameters: its amplitude D and the cutoff distance r c . 
If the molecular axis at a lattice site has three orientations, say the x, y or z axes, then when D = 0, 
the system is equivalent to the 3-state Potts model: the transition to the disordered phase is known 
to be of first order. When the ground-state configuration is shown to be composed of two 

independent interpenetrating layered subsystems which form a sandwich whose periodicity depends 
on D and r c . We show by extensive Monte Carlo simulation with a histogram method that the 
phase transition remains of first order at relatively large values of r c . 

PACS numbers: 61.30.Cz, 64.70. M-: 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

The nature of the transition from one phase to an- 
other is one of the most important problems in statisti- 
cal physics and in various areas of materials science not 
limited to physics. Since the introduction of the renor- 
malization groupAiS, with new concepts based on the sys- 
tem symmetry and scale, the understanding of the na- 
ture of the phase transition in many systems becomes 
clear. However, there are complicated systems where 
that method encounters much of difficulties in applica- 
tion. One can mention frustrated systems where com- 
peting interactions cause highly degenerate ground states 
with instabilities^. 

In this paper, we are interested in the phase transi- 
tion in molecular crystals which has been a subject of 
intensive investigations since 40 years. This spectacular 
development was due to numerous applications of liquid 
crystals in daily life^£. Liquid crystals are somewhere be- 
tween solid and liquid states where molecules have some 
spatial orientations which, under some conditions, can 
order themselves into some structures such as nematic 
and smectic phases. Liquid crystal display uses prop- 
erties of these ordered phases^. In spite of the large 
number of applications using experimental findings, the- 
oretical understanding in many points is still desirable. 
One of the most important aspects of the problem is the 
origin of layered structures observed in smectic order- 
ing: in smectic phases the molecules are ordered layer 
by layer with periodicity. Experiments have discovered, 
for example, smectic phases of 3-layer, 4-layer, or 6-layer 
periodicity— Different theories have been suggested to 
interpret these observations^—. One of the unanswered 
questions is what is the origin of the long-period lay- 



ered structure in observed smectic phases? Hamaneh et 
ali^i suggested that the effective long-range interaction is 
due to bend fluctuations of the smectic layers which may 
stabilize commensurate structures in particular the six- 
layer phase. This suggestion is too qualitative to allow a 
clear understanding of the long-period layered structure 
in smectic phases. To our opinion, we should distinguish 
long-range correlation and long-range interaction. We 
know that long-range correlation can be a consequence 
of short-range interactions, not necessarily of long-range 
ones, as it is well known in theory of critical phenomena 2 . 
However, a long-range correlation does not imply a par- 
ticular long-period configuration. 

To look for a physical origin of long-period structures, 
we concentrate ourselves in the present paper to the ef- 
fect of a dipolar interaction in a Potts model, in addition 
to an exchange interaction between nearest neighbors 
(NN). We suppose that each molecule has a molecular 
axis which can lie on one of the three principal direc- 
tions. An example of such a molecule is the ammoniac 
molecule NH3. Without the dipolar interaction, the in- 
teraction between neighboring molecules gives rise to an 
orientational order of molecular axes at low temperature. 
In this situation, the system can be described by a 3-state 
Potts model in three dimensions (3D). The Potts model 
is very important in statistical physics. It is in fact a class 
of models each of which is defined by the number of states 
q that an individual particle or molecule can have. The 
interaction between two neighboring molecules is nega- 
tive if they are in the same state, it is zero otherwise. 
Exact solutions are found for many Potts models in two 
dimensions (2D)22. In 2D, the phase transition is of sec- 
ond order for q < 4 and of first order for higher q. In 3D, 
a number of points are well understood. For example, the 
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phase transition is of first order for q > 2. The molecular 
crystal described above undergoes therefore a first-order 
transition in the absence of a dipolar interaction. 

The purpose of this paper is to investigate the effect of 
the dipolar interaction on the ground state (GS) struc- 
ture and on the nature of the phase transition. The dipo- 
lar interaction is a long-range interaction which yields 
different GS structures depending on the shape of the 
sample as will be discussed in the next section. To carry 
out our purpose, we use a steepest descent method for the 
GS determination and the Monte Carlo (MC) simulation 
combined with the histogram technique to distinguish 
first- and second-order characters. 

In section |nl we show our model and analyze the GS. 
Results of MC simulations are shown and discussed in 
section Hill Concluding remarks are given in section HV1 

II. MODEL AND GROUND-STATE ANALYSIS 

We consider a simple cubic lattice where each site is 
occupied by an axial molecule. The molecular axis can be 
in the x, y or z direction. Let us denote the orientation 
of the molecule at the lattice site i by a unit segment, 
not a vector, which can lie in the x, y or z direction. We 
attribute the Potts variable cr=l, 2 or 3 when it lies in 
the x, y or z axes, respectively, in the calculation. By 
the very nature of the model shown below, the results do 
not depend on these numbers as in the standard g-state 
Potts model. 

We suppose that the interaction energy between two 
nearest molecules is - J (J > 0) if their axes are parallel, 
zero otherwise. With this hypothesis, the Hamiltonian is 
given by the following 3-state Potts model: 

H = -J^% ilffj ) (1) 

(id) 

where Oj is the 3-state Potts variable at the lattice site i 
and i s made over the nearest sites Oi and Uj. 

The dipolar interaction between Potts variables is writ- 
ten as 

Hd = D ~ 3 [S(<Ii) ■ U -f^ ] ■ U - ] Wat, a,) 

(id) r « ^ 

(2) 

where Ujj is the vector of unit length connecting sites 
i and j, D a positive constant depending on the mate- 
rial, S(c7j) is defined as the unit vector lying on the axis 
corresponding to the value of <ii. The sum is limited at 
some cutoff distance r c . Without the Kronecker condi- 
tion and in the case of classical XY or Heisenberg spins, 
the dipolar interaction gives rise to spin configurations 
which depend on the sample shape. For example, in 2D 
or in rectangular slabs spins lie in the plane to minimize 
the system energy. 

Let us first discuss about the GS in the Potts model in- 
troduced above. When D = the GS is uniform with one 



orientation value, namely it is 3-fold degenerate. How- 
ever, for a nonzero D, the ground state changes with 
varying r c . To determine the GS, we use the steepest 
descent method which works very well in systems with 
uniformly distributed interactions. This method is very 
simple 23,24 (i) we generate an initial configuration at ran- 
dom (ii) we calculate the local field created at a site by 
its neighbors using (fTJ) and @ (iii) we take the Potts 
variable to minimize its energy (i. e. align the "Potts 
spin" in its local field) (iv) we go to another site and re- 
peat until all sites are visited: we say we make one sweep 
(v) we do a large number of sweeps per site until a good 
convergence is reached. 

We show some examples in Fig. Q] For very small D, 
the GS is uniform for any r c as shown in Fig. QJ,. For 
increasing D 7 the GS has layered structures with period 
p = 1 (alternate single layers), p — 2 (double layers), 
p = 3 (triple layers),... depending on r c and D. 

We note that in Fig. [T] there are only two kinds of 
molecular orientations represented by two colors (on line) 
in spite of the fact that we have three possible orien- 
tations. Let us discuss about the single-layer structure 
shown in Fig. QJ. It is very important to note that 
each layer has no coupling with two nearest layers whose 
molecules lie on another axis. However, it is coupled to 
two next-nearest layers of the same color, namely the 
same molecular axis, lying within r c . In other words, the 
GS is composed of two interpenetrating "independent" 
subsystems. The same is true for other layered struc- 
tures: in a double-layer structure molecules of the same 
orientation (same color on line) interact with each other 
but they are separated by a double layer of molecules 
of another orientation. To our knowledge, this kind of 
GS has never been found before. It may have important 
applications at macroscopic levels. 

Note also that if the molecules of the single-layer struc- 
ture are successively in x-oriented and y-oriented planes, 
then the stacking direction of these planes is the z direc- 
tion. The molecules can choose the x and z orientations 
or the y and z directions. In those cases the stacking 
directions are respectively the y and x directions. These 
three possibilities are equivalent if the sample is cubic. 
In rectangular shapes, the stacking direction is along the 
smallest thickness. 

In order to understand how such GS configura- 
tions found by the steepest descent method depend 
on D and r c , we have considered the structures 
121212. ..(single-layer structure), 11221122... (2-layer 
structure), 111222111222... (3-layer structure) and car- 
ried out the calculations of the energy of a molecule Oi 
interacting with its neighbors Uj. The case where the 
configuration is uniform, i. e. there is only one kind of 
molecular orientation, say axis x, the dipolar energy of 
o-j is 

1 (v x -) 2 

^e^-3^h (3) 

Note that with the use of the Potts model for the dipolar 
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term in Eq. ((2J, the energy depends only on the axis, 
not on its direction as seen by the square term (h?-) 2 . 
In this sense, the model is suitable to describe axial, but 
non-directed, interacting molecules. 

If we transform the sum in integral, the sum in the 
first term gives 47rlnr c (integrating from 1 to r c ), while 
the second term gives -47r In r c , which cancels the first 
term. This is valid for r c larger than 1. Thus the dipolar 
energy Ed — D [4w In r c — Air In r c ] is zero for the uniform 
configuration in 3D space in the framework of our model. 
The energy of the system comes from the short-range 
exchange term, Eq. (fTJ). 

The energy per site of other layered structures with pe- 
riodicity p = 2 and 3 is numerically calculated and shown 
in Fig. E| for D/J = 0.8 and 2. In each figure, the GS 
is the structure which corresponds to the lowest energy. 
One sees in Fig. [2^ the following GS configurations with 
varying r c : 

- uniform GS: for 1 < r c < 1.3 (zone I) 

- single-layer structure: for 1.3 < r c < 1.8 (zone II) 

- double-layer structure: for 1.8 < r c < 3.65 (zone III) 

- triple-layer structure: for r c > 3.65 (zone IV) 

For D/ J — 2 there are only two possible GS with vary- 
ing r c as shown in Fig. [2Jd. 

We summarize in Fig. [3] the different GS in the space 
(r c , D) with J = 1. We note that for a given D, for ex- 
ample D — 0.8, the GS configuration starts with uniform 
configuration then with period 1, 2, 3, ... for increasing 
r c . At large D, uniform configuration is not possible at 
any r c . Long-period configurations are on the other hand 
favored at small D and large r c . 



III. PHASE TRANSITION: RESULTS 

We consider a sample size of N x N x N z where N and 
N z vary from 24 to 48 but N z can be different from N 
in order to detect the shape-dependence of the ground 
state.. The 3-state Potts model with NN exchange in- 
teraction J = 1 is used to describe the three molecular 
orientations. For the dipolar term, a cutoff distance r c is 
taken up to VTTj — 3.16 lattice distance. At this value, 
each molecule has a dipolar interaction with 146 neigh- 
bors. Periodic boundary conditions in all directions are 
employed. 

We have used the standard MC method 25 with the 
system size from 24 3 to 48 3 . The equilibrating time N\ 
is about 10 6 MC steps per site, and the averaging time N% 
is between 10 6 and 10 7 MC steps per site. The averages 
of the internal energy < U > and the specific heat Cy 
are defined by 



(u) = <n- 
(u 2 ) 



c 



v 



H d > 

-(U) 2 



(4) 
(5) 




k B T 2 

where < ... > indicates the thermal average taken over 
N2 microscopic states at T. 



FIG. 1. Ground state in the case where a) D/J — 0.4, r c = 
2.3: uniform configuration ; b) D/J = 1, r c = 2.3: double- 
layer structure; c) D/J = 2, r c = 2.3: single-layer structure; 
d) D/J — 0.4, r c = \/T0 ~ 3.16: triple-layer structure. See 
text for comments. 



We define the Potts order parameter Q by 

Q = [ ? max(Q 1 , Q 2 , Q q ) - l]/(q - 1) (6) 
where Q n is the spatial average defined by 

Q n - ~ xNxN *) (7) 
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FIG. 2. Ground state energy versus r c in the case where (a) 
D/J = 0.8 (b) D/J — 2. Blue, black, red and green lines 
represent GS energy of the uniform, single-layer, double-layer 
and triple-layer structures, respectively. Zones I, II, II, IV 
indicate these respective different GS. See text for comments. 
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FIG. 3. Ground state in space (r c ,D): the numbers 0, 1, 
2 and 3 denote the uniform, single-layer, double-layer and 
triple-layer structures, respectively. The first column displays 
the values of r c with the number of neighbors indicated in the 
parentheses. 



n(n = l,.,.,q) being the value of the Potts variable at 
the site j. For q = 3, one has n — 1,2,3 representing re- 
spectively the molecular axis in the x, y and z directions. 
The susceptibility is defined by 



X = 



(Q 2 ) - {Q? 



(8) 



In MC simulations, we work at finite sizes, so for each 
size we have to determine the " pseudo" transition which 
corresponds in general to the maximum of the specific 
heat or of the susceptibility. The maxima of these quan- 
tities need not to be at the same temperature. Only 
at the infinite size, they should coincide. The theory of 
finite-size scaling 2 ^— permits to deduce properties of a 
system at its thermodynamic limit. We have used in this 
work a size large enough to be close to the extrapolated 



bulk transition temperature. 

We show first in Fig. [4] the energy per site E =< U > 
/(N x N x N z ) and the order parameter M =< Q > 
versus T, for several lattice sizes in the absence of the 
dipolar interaction, i. e. D = 0. As said earlier this case 
corresponds to the 3-state Potts model. We find a very 
sharp transition. In order to check the first-order na- 
ture of this transition, we used the histogram technique 
which is very efficient in detecting weak first-order tran- 
sitions and in calculating critical exponents of second- 
order transition o 27 ' 28 . The main idea of this technique 
is to make an energy histogram at a temperature To as 
close as possible to the transition temperature. Often, 
one has to try at several temperatures in the transition 
region. Using this histogram in the formulae of statis- 
tical physics for canonical distribution, one obtains en- 
ergy histograms in a range of temperature around To. In 
second-order transitions, these histograms are gaussian. 
They allows us to calculate averages of physical quantities 
as well as critical exponents using the finite-size scaling. 
In first-order transitions, the energy histogram shows a 
double-peak structure at large enough lattice sizes. Us- 
ing this method, we have calculated the histogram shown 
in Fig. [5] for several lattice sizes. With increasing size, 
the two peaks are well separated with the dip going down 
to zero, indicating a tendency toward an energy discon- 
tinuity at the infinite size. The distance between the two 
peaks is the latent heat AT. 




FIG. 4. Internal energy per spin E (upper) and order parame- 
ter M =< Q > (lower) versus temperature T without dipolar 
interaction, for several sizes N z — N — 24 (black circles), 36 
(void circles), 48 (blue triangles, color on line), J = 1. 

We take into account now the dipolar interaction. In 
order to see its progressive effect, we perform simula- 
tions with r c varying between r c = y/6 and r c = ^/lQ 
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FIG. 5. Energy histogram at the transition temperature with- 
out (D = 0, upper) and with (D — 2, lower) dipolar interac- 
tion, for several sizes N z — N = 24 (black circles), 36 (void 
circles), 48 (blue triangles, color on line), r c — \/T0, J = 1. 



and follow the change of the characteristics of the phase 
transition. 

We show in Fig. [6]the energy and the order parameter 
versus T for D = 2, r c = vTU. The transition is still very 
sharp but the strong first-order character diminishes. 

We show in Fig. [7] the energy histogram for several 
values of r c with D = 2. As seen the latent heat becomes 
small at r c = y/lO. 

We observe that the latent heat AE diminishes with 
increasing r c . However, we cannot conclude that the 
first-order disappears at large r c . To check that point 
we need to go to larger r c which will take a huge CPU 
time because of the increasing number of neighbors (we 
recall that for r c = vlO, we have 146 neighbors for each 
molecules) . We observe from Fig. [7] that the latent heat 
does not change significantly with r c up to r c ~ 2.85, 
meaning that the first order transition is dominated by 
the short-range 3-state Potts interaction. The latent heat 
decreases rather strongly afterward but we do not know 
if it tends to zero or not. 



IV. CONCLUSION 

We have studied in this paper a molecular crystal char- 
acterized by three possible orientations of the molecular 
axes. The model is described by a short-range 3-state 
Potts model and a dipolar Potts interaction. We have 
analyzed the ground state as functions of the dipolar 
interaction strength D and its cutoff distance r c . The 
ground state shows remarkable properties: it is composed 
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FIG. 6. Internal energy per spin E and order parameter 
M =< Q > versus T for D — 2, r c = \/T0 with several 
lattice sizes N z = N — 24 (black circles), 36 (void circles), 48 
(blue triangles, color on line), and J = 1. 
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FIG. 7. Energy histogram showing double-peak structure for 
D — 2 at several values of r c : from left to right r c = v6 
(black circles), y8 (void circles), 3 (red stars, color on line) 
and vlO (blue triangles, color on line). N z — N = 48, J = 1. 



of two interpenetrating independent subsystems with a 
layered structure whose periodicity depends on D and 
r c . To our knowledge, such a ground state has not been 
seen in uniformly interacting molecules. The stacking 
of independent ordered layers is reminiscent of a smec- 
tic ordering. Note that the model is applied to smectic 
phases where tilted angles are uniform: it can be of type 
A as supposed here or of type C if the molecular axis 
does not coincide with the x, y or z axis of the lattice. 
The coupling between adjacent different ordered planes 
is here strictly zero at T = 0, due to the Potts condition 
in the model. We have used the Monte Carlo histogram 
method to study the phase transition in this system. In 
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the absence of the dipolar interaction, the model which 
is a 3-state Potts model is known to undergo a first-order 
phase transition from the orientational ordered phase to 
the disordered state. Upon the introduction of a dipolar 



interaction of sufficient strength into the Potts model, the 
ground state is broken into layers as described above. We 
have shown that the transition remains of first order as 
the cutoff distance is increased at least up to r c = \/l0. 
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